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摘要 :【 目的 】 本 研究 旨 在 筛选 小 菜 蛾 Plutella xylostella s] 3 98 EF RIE Isaria fumosorosea 侵 染 的 
免疫 应 答 及 其 网 络 调控 基因 ,以 进一步 探讨 小 菜 蛾 对 玫 烟 色 棒 束 孢 的 防御 机 制 。【 方 法 ] 采 用 第 二 
代 高 通 量 测序 技术 RNA-seq, 对 处 理 后 12 h 的 感染 玫 烟 色 棒 束 孢 和 健康 (平行 对 照 ) 的 小 菜 蛾 4 龄 
幼虫 转录 组 进行 了 测序 ,通过 生物 信息 学 分 析 对 得 到 的 差异 基因 进行 了 功能 注释 和 分 类 ,对 差异 基 
因 参 与 的 信号 通路 进行 了 分 析 。【 结果】 在 感 菌 和 健康 小 菜 蛾 幼虫 转录 组 两 个 表达 谱 里 ,分 别 获 得 
了 12 346 987 和 12 315 210 个 clean reads, 有 60.93% fe 61. 2696 &*] reads 分 别 能 比 对 到 参考 基因 库 
里 ,其 中 完美 匹配 (perfect match) 的 比例 分 别 为 32. 1596 和 32.73% 。 共 得 到 351 个 显著 差异 表达 
基因 (differentially expressed unigenes, DEUs) ,上 调 表达 基因 有 275 个 ,下 调 表达 基因 有 76 个 ,与 免 
疫 防 御 反 应 潜在 相关 的 基因 有 156 个 。GO (Gene Ontology) 富 集 分 析 有 102 个 DEUs 分 布 到 46 个 
GO term 里 ,KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway 富 集 分 析 结 果 显 示 , 有 132 
个 DEUs 显著 富 集 在 13 个 代谢 通路 (pathway) 里 。 [2536] 这 些 差异 表达 基因 中 ,大 部 分 编码 潜在 
的 与 免疫 识别 及 调控 相关 的 基因 ,集中 在 能 量 代 谢 、 疾 病 反应 和 防御 反应 等 相关 通路 。 研 究 结果 为 
挖掘 与 玫 烟 色 棒 束 孢 侵 染 相关 的 小 菜 蛾 免疫 应 答 候选 基因 提供 了 重要 数据 库 , 也 为 益 明 小 菜 蛾 对 
玫 烟 色 棒 束 孢 的 免疫 机 制 黄 定理 论 基 础 。 
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Abstract: [Aim] This study aims to explore the immune-related genes and immune recognition 
regulatory network of Plutella xylostella against the entomopathogenic fungus /saria fumosorosea, so to 
further investigate the innate immune mechanisms of P. xylostella to I. fumosorosea. [ Methods] We 
analyzed the transcriptome of I. fumosorosea-infected and uninfected (the control) 4th instar larvae of P. 
xylostella at 12 h after treatment using high-throughput single-end RNA-sequencing ( RNA-seq). The 
differentially expressed genes and their functions, classifications and signaling pathways were analyzed 
using bioinformatic tools. [Results] The numbers of 12 346 987 and 12 315 210 clean reads were obtained 
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for the infected and uninfected libraries, respectively. Among the two RNA-seq libraries, 60. 93% and 
61.2696 of reads could be mapped to reference genome with a perfect match proportion of 32. 1596. and 
32.73% , respectively. A total of 351 differentially expressed unigenes ( DEUs) were identified between 
infected and uninfected larvae. Among them, 275 DEUs were up-regulated while 76 DEUs down- 
regulated. In total, there were 156 DEUSs putatively involved in immune responses. A Blast2GO analysis 
indicated that 102 differentially expressed transcripts were able to be assigned to 46 Gene Ontology ( GO) 
terms. Kyoto Encyclopedia of Genes and Genomes ( KEGG) pathway enrichment analysis mapped 132 
DEUs to 13 pathways ( P «0. 05). [Conclusion] Most of the DEUs are putatively linked to innate immune 
recognition and regulation. These infection-responsive DEUs are enriched in various immune processes, 
such as energy metabolism, disease and defense response. The findings generated from this study provide 
a bioinformatic database for searching immune-related genes, and also benefit the future investigations of 
the innate immune mechanisms of P. xylostella against the entomopathogenic fungus 1. fumosorosea. 
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ANS Plutella xylostella 是 一 种 为 害 十 字 花 科 
蔬菜 的 世界 性 农业 害虫 (Talekar and Shelton, 
1993) , 目前 对 用 于 防治 其 危害 的 各 类 化 学 杀 虫 剂 
都 能 较 快 产生 抗 性 (Mohan and Gujar, 2003; Baek et 
al., 2005) ,其 亦 是 最 早 被 发 现在 田间 对 Bt 制剂 产 
生 抗 性 的 害虫 ( Heckel et al., 1999) ,人 研究 新 的 防治 
途径 ,利用 昆虫 病原 真菌 进行 生物 防治 是 持久 高 效 
控制 小 菜 蛾 为 害 的 必然 趋势 ( BAR SC, 2015; 雷 仲 
仁 等 , 2016) , 

SUI GPS YR Yl. [saria fumosorosea 作为 一 类 重要 
的 生 防 真菌 ,能 寄生 多 种 蔬菜 和 果树 害虫 ( D 00 d 
TIO HH J6, 1999; Wraight et al., 2000; Meyer et al., 
2008; Marjorie et al., 2010) , X 4E EK f& À 1E 
Beauveria bassiana ÑM & f T 218 W Metarhizium 



































与 抗 免 疫 的 动态 互 作 关系 ( 雷 妍 圆 等 , 2011), 
宿主 免疫 系统 与 昆虫 病原 真菌 的 作用 是 多 方面 
的 ,昆虫 的 生理 防御 基因 在 不 同 微生物 诱导 下 的 表 
达 模 式 不 同 ( Wang et al., 2010) ,昆虫 病原 真 苗 进入 
虫 体 之 后 ,被 昆虫 的 天 然 免疫 受 体 即 模式 识别 受 体 
( pattern recognition receptors, PRRs) 识别 ,开始 免疫 
反应 。 虫 体 天 然 免 疫 系 统 对 病原 菌 的 抵抗 ,是 延缓 
病原 真菌 致 病 速 度 的 主要 原因 。 有 研究 表明 ,人 为 
抑制 昆虫 的 免疫 系统 , 其 生存 将 受到 影响 
( Magalhaes et al., 2010 ) 。 找 到 免疫 应 答 过 程 中 的 
关键 因子 并 以 此 为 靶 标 , 筛选 对 这 个 关键 因子 有 抑 
制 \ 阻 断 等 功能 的 基因 ,对 于 高 效 杀 虫 微生物 的 篇 
选 改良 和 设计 具有 十 分 重要 的 应 用 价值 。 
目前 小 汪 蛾 对 玫 烟 色 棒 束 孢 免疫 受 体 基因 及 免 










































































anisopliae 之 后 最 有 开发 潜力 的 昆虫 病原 真 
( Milner, 1997; Wraight et al., 2000) , F] FH ECH £5, 
棒 束 孢 控制 虫害 已 得 到 广泛 认同 ,但 目前 国内 仍 未 
实现 该 菌 的 产业 化 生产 和 田间 应 用 ,除了 其 发 酵 生 
产 工艺 和 施用 制剂 还 有 待 研 发 ,来 自 寄主 昆虫 体内 


























疫 调控 网 络 关键 基因 的 种 类 尚 不 清楚 ,这 些 免疫 相 
关 基 因 的 抗菌 功能 和 互 作 关系 网 有 待 进一步 的 研 
究 。 为 此 ,本 研究 从 生物 信息 学 人 手 , 利 用 RNA-seq 
表达 谱 技 术 ( Wang et al., 2009) ,对 感 菌 与 健康 小 菜 
我 进行 表达 谱 测 序 和 差异 分 析 , 比较 感 落 与 健康 小 菜 
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的 抵抗 作用 也 成 为 了 影响 其 防 效 的 一 个 重要 因素 。 
本 研究 室 曾 报道 一 株 高 致 病 力 玫 烟 色 棒 束 孢 ( A 
华 等 , 2007) , E 1x 10 孢子 /mL 浓度 下 对 小 菜 蛾 2 
龄 幼虫 的 LTs 仅 1.72 d。 我 们 对 该 菌株 与 小 菜 蛾 的 
互 作 关系 进行 了 较为 系统 的 研究 ,发 现 小 菜 蛾 4 龄 
幼虫 感 菌 24 h 后 , 虫 体 血细胞 总 数 由 正常 的 每 mm 
血 淋 巴 中 750 个 大 幅 上 升 至 最 高 值 6 250 个 , 主要 
以 具有 人 免疫 功能 的 浆 血 细胞 和 粒 血细胞 对 菌 丝 段 进 
行 否 哈 千 包 。 面 对 虫 体 的 免疫 识别 与 防御 , 菌 体 
采取 了 增殖 小 型 颗粒 状 菌 丝 段 快速 侵 染 的 逃避 策 
略 , 昆 虫 与 病原 真菌 之 间 存 在 着 识别 与 反 识别 免疫 























蛾 差异 基因 表达 量 的 变化 ,寻找 与 玫 烟 色 棱 束 孢 侵 染 
相关 的 小 菜 蛾 免疫 应 管 候选 基因 ,为 进一步 探讨 特异 
性 表达 差异 基因 与 免疫 的 关系 奠定 信息 学 基础 。 











1 材料 与 方法 


1.1 供 试 昆虫 和 真菌 

本 试验 所 用 的 小 菜 蛾 来 自 于 华南 农业 大 学 农学 
院 天 敌 昆 虫 实验 室 , 为 人 工 饲 养 多 代 的 室内 种 群 , 室 
内 温度 控制 为 (25 +1) ,相对 湿度 60% ~90% , J6 
周期 14L: 10D ,幼虫 长 至 4 龄 初 供 试 。 政 烟 色 棱 束 
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fü Tk EBCL 03011 来 源 于 台湾 乳白 蚁 Coptotermes 
formosanus 染 菌 尸体 ( 雷 妍 圆 等 , 2011) , 菌 种 于 华南 
农业 大 学 昆虫 生态 室 -800 下 保藏 。 菌 株 在 萨 氏 培 
养 基 SDAY( 葡萄糖 40 g, EAK 10 g ,酵母 膏 10 g, 
琼脂 15 ~20 g, Zik 1 L) 平板 上 ,25% ,无 光照 条 
件 下 培养 20 d ,用 接种 环 将 分 生 孢 子 粉 刮 至 灭 菌 的 
0.0596 Tween-80 溶液 中 ,配制 成 浓度 为 1 x 10° $8 
Tal 的 孢子 悬浮 液 供 试 。 
1.2 小 菜 蛾 的 免疫 诱导 及 表达 谱 测 序 样品 制备 

用 软 毛 笔 挑 选 个 体 大 小 活力 一 致 的 4 龄 幼虫 ， 
分 感 菌 和 健康 对 照 两 部 分 处 理 。 一 部 分 用 孢子 悬浮 
液 浸渍 10 s 后 挑 出 , 另 一 部 分 用 灭 菌 的 0.05% 吐 
温 -80 溶液 浸渍 10 s 处 理 ,分 别 将 各 处 理 的 幼虫 置 
于 滤纸 上 , 吸 去 多 余 水 分 后 移 至 新 鲜 甘蓝 叶片 上 , 放 
AERA HE TE ACT EE LC 直径 =9 cm) 中 ， 
培养 下 用 保鲜 膜 封 口 保湿 ,封口 处 用 针 扎 数 孔 以 通 
风 透 气 。 幼 虫 接种 后 放 在 250C , 光 周 期 14L: 10D 的 
人 工 智能 箱 中 饲养 。 收 集 感 菌 12 h 的 小 菜 蛾 幼虫 
作为 感 菌 样品 ,平行 收集 0.05% 吐 温 -80 溶液 处 理 
12 h 的 健康 幼虫 为 对 照样 品 ( 雷 妍 圆 等 , 2011) 。 将 
存 有 样品 的 离心 管 置 于 液 氮 中 速冻 ,然后 保存 于 
-80'C Bj IE RNA 降解 。 每 个 样品 收集 7 ~8 头 幼 
E ,提取 RNA 用 于 测序 。 
1.3 RNA 抽 提 、cDNA 文库 构建 及 solexa 测序 

按 Trizol 试剂 抽 提 法 提取 总 RNA, DA 1% 琼脂 
糖 凝 胶 电 泳 检 测 RNA 的 完整 性 ,用 Agilent 2100 检 
i| RNA 的 纯度 。 用 带 有 0ligo (dT) 的 磁 珠 富 集 含 
有 poly-A 的 mRNA ,然后 通过 Fragmentation buffer 
对 富 集 到 的 mRNA 进行 片段 化 处 理 ,以 随机 打 断 成 
短片 段 的 mRNA 作为 模板 ,再 通过 六 碱 基 的 随机 引 
物 (random hexamers) 合成 cDNA 第 1 链 ,随后 加 入 
DNA polymerase 工 缓冲 液 和 RNase H 合成 cDNA 第 
2 链 , 用 QiaQuick PCR. ORAE cDNA 双 链 进行 纯 
化 ,经 过 未 端 修复 、 加 碱 基 A、 加 测序 接头 等 一 系列 
步骤 ,用 于 脂 糖 凝 腕 电泳 回收 目的 片段 ,最 后 通过 
PCR 扩 增 目的 片段 完成 cDNA 文库 的 构建 。 合 成 的 
cDNA 样品 由 深圳 华 大 基因 研究 院 (BGI) 进行 测序 ， 
测序 平台 为 Humina HiSeq 2000, 每 个 样品 1 个 
Line 的 测序 量 。 
1.4 原始 数据 处 理 

测序 产生 的 原始 图 像 数 据 通过 pase calling 转 
换 为 序列 数据 (raw data)。 利 用 固定 脚本 对 raw 
data 进行 一 系列 的 前 处 理 ,过 滤 掉 只 含有 测序 接头 
(adapter) 的 reads, 含 有 空白 片段 N 的 比例 超过 



















































































































































































10% 的 reads, 以 及 质量 值 Q «5. 碱 基数 目 比 例 超过 
整个 read 5096 以 上 的 低 质量 reads ,获得 可 用 于 分 析 
的 数据 Clean reads。 同 时 ,应 用 Trinity 软件 进行 混 
合 拼接 ,将 两 个 测序 样本 中 的 clean reads 与 玫 烟 色 
束 孢 转录 组 ( 竺 发表) 进行 比 对 ,去 除 虫 体 中 含有 
的 玫 烟 色 棒 束 孢 带 来 的 物种 序列 污染 ,得 到 非 元 余 
的 Unigene。 
1.5 Unigene 与 参考 序列 的 比 对 

将 各 样本 拼接 得 到 的 Unigene 用 SOAP 软件 (Li 
et al., 2009) 比 对 到 小 菜 蛾 基因 组 上 (http://iae. 
fafu. edu. en/DBM/index. php) ,产生 注释 。 
1.6 显著 差异 表达 基因 筛选 

分 别 计 算 感 菌 和 健康 小 菜 蛾 两 个 样本 中 所 有 基 
因 的 表达 量 , 对 样本 间 的 基因 表达 水 平 进行 差异 比 
较 。 基 因 表 达 量 的 计算 方法 使 用 RPKM 法 (reads 
per kb per million reads) ( Mortazavi et al., 2008) ,其 
计算 公式 为 : 
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假设 目标 基因 A 的 表达 量 为 RPKM ( A) , Wu C 
为 唯一 比 对 到 该 基因 A 上 的 clean reads 的 数目 ,NN 
为 唯一 比 对 到 参考 基因 序列 上 的 总 的 clean reads 的 
数目 ,L 为 目标 基因 A 的 总 碱 基数 目 。 使 用 edgeR 
( Robinson et al., 2010) 进行 差异 分 析 。 通 过 控制 
FDR(false discovery rate) [HR kE P Hh BE, 
行 多 重 检 验 , FDR 不 能 超过 5% 。 因 此 ,以 FDR < 
0. 00171 | log; Ratio | (1. 作为 闪 值 判断 基因 表达 差异 
的 显著 性 ,筛选 出 显著 差异 表达 基因 (differentialty 
expressed unigenes, DEUs) 。 
1.7 GO 功能 注释 和 KEGG pathway 富 集 分 析 

将 DEUs 比 对 到 GO 数据 库 和 KEGG 数据 库 。 
利用 Blast2GO 对 DEUs 进行 GO 功能 注释 , 通 i 
WEGO( Ye et al., 2006) 对 DEUs 在 生物 学 过 程 、 
子 功能 和 细胞 组 分 3 个 部 分 的 功能 和 作用 方式 进行 
GO 分 类 。 通 过 MADIBA ( Law et al., 2008 ) 分 析 
DEUs 所 处 的 代谢 路 径 或 信号 传导 途径 。 运 用 超 几 
何 检验 的 方法 ,确定 GO M KEGG pathway 在 DEUs 
中 的 P 值 , 当 P<0.05, 则 认为 基因 在 特定 GO 或 
pathway 中 是 显 车 富 集 的 ,反之 则 不 富 集 。 


RPKM = 
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2 结果 


2.1. 测序 结果 分 析 
从 健康 和 感 菌 小 全 蛾 转录 组 中 分 别 获 得 了 
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12 565 677 和 12 599 166 个 raw reads (图 1) 。 除 去 
含有 接头 的 、 含 未 知 碱 基 和 低 质量 的 reads, 分 别 获 
f& f 12315 210 和 12 346 987 个 clean reads( 表 1)， 
对 所 得 到 的 clean reads, 使 用 reads 比 对 软件 
SOAPaligner/soap2 将 clean reads 分 别 比 对 到 小 菜 蛾 
基因 组 上 上。 其 中 ,达到 完美 匹配 的 reads 分 别 为 健 
康 4 031 042 条 , 感 菌 3 969 192 条 ,分 别 有 32. 7396 
和 32.15% 的 reads 完美 匹配 到 小 菜 蛾 基因 组 数据 
中 唯一 的 unigene, 各 处 理 样 品 有 约 5096 的 reads 能 
唯一 匹配 到 小 荣 蛾 基因 组 数据 中 的 unigene( 允许 有 
2 bp 的 错 配 ) 。 


表 1 健康 和 感染 玫 烟 色 棒 束 孢 的 小 菜 蛾 RNA-seq 测序 结果 
Table 1 Statistics of RNA sequencing metrics from 
transcriptomes of Plutella xylostella larvae uninfected 


and infected by Isaria fumosorosea 





























度量 标准 健康 小 菜 蛾 感 菌 小 菜 蛾 
Metric Uninfected Infected 
Reads 总 数 Total reads 12 315 210 12 346 987 
碱 基 对 总 数 Total base pairs 603 445 290 605 002 363 
符合 设 定 参 数 的 reads 总 数 7 544 209 7 523 304 
Total mapped reads (61.26% ) (60.93% ) 
4 031 042 3 969 192 
r— 3 
完美 匹配 Perfect match (32.73% ) (32.15%) 
6 261 056 6 236 958 
L i h 
EGRE Umigoe wale (50. 84% ) (50.51% ) 


A 





0.46%, 0.03% 


1.50% 
& &N Containing N (3 352) 
回 低 质量 Low quality (188 687) 


m 90%-100% (2 383) 
»" 80%-90% (1375) 
25 7096-8095 (1 363) 
心 60%-70% (1 332) 
1360) 
1328) 


人 | 1 15095-6096 
| 11140%-50% 
] * .*30%-40% (1 380) 
$$ 20%-30% (1 404) 
HH 10%-20% (1 460) 
£2: 0-10% (1 023) 


po NEHME S 





E 含 接头 Containing adaptor (58 428) 


国 过 小 后 Clean reads (12 315 21g 





2.2 ”测序 质量 评估 及 基因 覆盖 度 统计 

对 两 个 样本 测序 出 来 的 raw reads 进行 数据 处 理 
(图 1) ,在 健康 和 感 菌 小 菜 蛾 样本 中 , 只 含有 接头 的 
序列 分 别 占 raw reads 的 0. 4696 fl 0. 5196 ,没有 测 
出 的 序列 在 两 个 库 中 各 占 0.03% ,质量 低 的 测序 分 
别 占 1.50% 和 1.47% , 剩 下 可 用 于 后 续 分 析 数 据 的 
clean reads 分 别 占 两 个 库 raw reads 的 98. 01% 和 
98.00% 。 基 因 测序 履 盖 度 指 每 个 基因 被 reads fs 
的 百分比 ,其 值 等 于 在 基因 中 唯一 能 比 对 到 参考 基 
因 组 上 的 reads 的 碱 基数 目 与 该 参考 基因 的 所 有 碱 
基数 目的 比值 (unique mapping reads) 。 在 健康 和 感 
菌 小 菜 蛾 两 个 样本 中 ,基因 和 覆盖 度 在 90% ~ 10096 
之 间 的 reads 在 每 个 样本 中 是 最 多 的 ,分 别 有 2 383 
和 2 518 条 序列 ;在 80% ~90% 之 间 的 分 别 有 1 375 
和 1 535 条 ;覆盖 度 在 0 ~10% 之 间 的 reads 最 少 , 分 
别 为 1023 和 972 条 (图 1)。 
2.3 ”显著 差异 表达 基因 筛选 

从 图 2 可 以 看 出 ,健康 和 感 菌 小 菜 蛾 之 间 共 得 
到 15 089 个 差异 表达 基因 ,其 中 显著 差异 表达 351 
个 (FDR<0.001 and |log,Ratio| 三 1) ,上调 表 达 基 
275 个 ,下 调 表 达 基 因 76 个 ,上 调 DEUs 数量 是 
38893 倍 ( 图 2: A) ,与 免疫 相关 的 基因 有 156 
个 。 在 这 些 差异 表达 基因 中 ,表达 量 差异 倍数 达到 





0.5194. 0.0396 47% 









W 含 接头 Containing adaptor (64 196) 
& &N Containing N (3 325) 

回 低 质量 Low quality (184 658) 

国 过 滤 后 Clean reads (12 3456 987) 


m 90965-1009 (2 518) 
.* 8096-9096 (1 535) 
4 7096-8096 (1 381) 
"^ 60%-70% (1 392) 
"T 50%-60% (1342) 
1! 4096-5066 (1343) 
'.' 3096-4096 (1 382) 
M 2096-3096 (1399) 
+ 1096-2096 (1380) 
J 0-10% (972) 


图 1 健康 (A，C) 和 感染 玫 烟 色 棒 束 孢 的 (B,，D) 小 菜 蛾 幼虫 转录 本 基本 参数 


Fig. 1 


Basic quantitative transcriptome parameters of transcriptomes of Plutella xylostella larvae uninfected ( A, C) 


and infected (B, D) by /saria fumosorosea 
A, B; 测序 质量 评估 图 Classification of raw reads; C, D; AA s EE Aht Distribution of gene coverage. 
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10 倍 以 上 的 高 表达 DEUs 有 2 条 ,为 CCG010114. 1 间 有 73 条 ,全 部 为 上 调 表 达 基 因 , 涉 及 “ 抗 
(log,Ratio = 10. 84) 和 CCG007755. 1 ( log, Ratio = ( pxCecropin a Tae 

10.83) (El 2: B)。 大 部 分 的 基因 ( ~ 78. 6% ) 共 素 ”( gloverin) “AKARE (catalase) “细胞 色素 
276 条 表达 量 都 分 布 在 1 ~2 倍 之 间 , 有 200 条 上 调 P450” ls P450) 等 可 能 与 免疫 防御 反应 相 
表达 ,7 
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Number of differentially expressed unigenes 








6 条 下 调 表达 ;表达 量 差异 倍数 在 2~10 售 关 的 注释 ,其 中 有 12 个 DEUs 没有 注释 。 
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图 2 健康 和 感染 玫 烟 色 棒 束 孢 的 小 菜 蛾 幼虫 转录 组 中 差异 表达 基因 


Fig. 2 Differentially expressed unigenes ( DEUs) between transcriptomes of Plutella xylostella larvae uninfected and 














infected by /saria fumosorosea 
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2.4 差异 基因 GO 富 集 分 析 能 集中 在 基因 生物 学 过 程 本 体 的 “细胞 过 
对 差异 表达 基因 进行 GO 功能 注释 分 析 。GO ( cellular process) ,' gi fer 过 程 ”( metabolic process) , 


二 级 分 类 表明 ,有 102 个 DEUs 注释 到 46 个 GO 
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表达 基因 统计 Statistics of DEUs; B; 差异 表达 基因 表达 量 差异 倍数 分 布 Distribution of folds of change of expression levels of DEUs. 
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“应 激 反应 ”(response to stimulus) 和 ”单一 有 机 体 过 


term 中 (图 3) ,包括 22 个 生物 学 过 程 ( biological 程 ”( single-organism process ) , 细胞 组 分 的 “细胞 ” 


4A” 


processes) ,12 个 细胞 组 分 (cellular component ) fil 10 (cell) 和 ”细胞 部 件 ”(cell parts) ,分 子 功能 的 “结合 
个 分 子 功能 (molecular function) 。 大 部 分 的 基因 功 (binding) 和 ”催化 作用 ”(catalytic) 等 GO term 
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3 ”健康 和 感染 玫 烟 色 棒 束 孢 的 小 菜 蛾 幼虫 转录 组 中 差异 表达 基因 GO 二 级 分 类 


Fig. 3 Gene Ontology classification (level-2) of differentially expressed unigenes ( DEUs) between transcriptomes of 





Plutella xylostella larvae uninfected and infected by /saría fumosorosea 
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结合 GO 数据 库 , 将 DEUs 比 对 上 的 3 个 功能 注 
释 本 体 进行 三 级 分 类 ,得 到 若干 个 子 集 ( 表 2) 。 对 
感 菌 和 健康 小 菜 蛾 间 DEUs 可 能 参与 的 基因 功能 进 
行 聚 类 统计 分 析 , 发 现 unigene 显著 富 集 在 细胞 组 
分 亚 类 的 “ 胞 外 区 域 ”(extracellular region ) ( GO: 
0005576) 以 及 分 子 功能 亚 类 的 “氧化 还 原 酶 活性 ” 
( oxidoreductase activity) ( GO; 0016491) ) “细菌 细胞 
表面 结合 ”( bacterial cell surface binding) ( GO; 
0051635) “细胞 表面 结合 ”(cell surface binding) 
(GO: 0043498) 等 。 在 生物 学 过 程 亚 类 中 , 差异 表 
达 基 因 和 集中 在 :“ 应 激 反 应 ”( response to stress ) 
(GO; 0006950 ) “防御 反应 ”(defense response ) 
(GO: 0006952 ) “应 激 反 应 中 对 翻译 的 调控 ” 
(regulation of translation in response to stress) ( GO: 
0043555) , “应 激 反 应 中 对 翻译 起 始 的 调控 ” 
(regulation of translational initiation in response to 
stress) ( GO: 0043558 ) “细胞 自我 吞噬 调控 ” 
(regulation of autophagy ) ( GO; 0010506) ) “天 然 免 


1i 














疫 反 应 调控 ”( regulation of innate immune response ) 
(GO; 0045088 ) “免疫 系统 过 程 ”(immune system 
process) ( GO; 0002376 ) “天 然 免疫 反应 ”(innate 
immune response ) ( GO; 0045087 ) 和 “免疫 反应 ” 
(immune response ) ( GO; 0006955 ) 等 一 系列 与 免疫 
反应 相关 的 功能 ,说 明 小 菜 蛾 被 玫 烟 色 棱 束 孢 侵 染 之 
后 ,体内 产生 了 一 系列 的 生理 生化 防御 反应 ,主要 涉 
及 能 量 代 谢 、 转 录 和 翻译 调控 以 及 防御 反应 等 方面 。 
2.5 差异 表达 基因 KEGG pathway 富 集 分 析 

将 健康 和 感 菌 小 菜 蛾 间 显 闭 差 异 表达 基因 
DEUs 映射 到 KEGG 路 径 库 里 富 集 分 析 ( 表 3) , 结 
果 显 示 , 一 共有 163 个 DEUs 比 对 到 153 个 代谢 通 
路 (pathway) 里 ,其 中 有 132 个 DEUs 显著 富 集 在 13 
条 代谢 通路 中 (Q < 0. 05 ) ,说 明 这 些 通路 在 真菌 诱 
导 下 被 显著 调控 。 富 集 DEU 最 多 的 通路 为 “A 型 
流感 “Influenza A) 代谢 通 路 (12 个 上 调和 8 个 下 
调 DEUs) “和 蛋白 消化 和 吸收 代谢 途径 ”(Pprotein 
digestion and absorption) (10 个 上 调和 8 个 下 调 DEUs) , 




























































































表 2 ， 健康 和 感染 玫 烟 色 棒 束 孢 的 小 菜 蛾 幼虫 转录 组 中 差异 表达 基因 ( 己 值 <1) GO 三 级 分 类 
Table2 Significantly enriched GO term (level-3) (P-value <1) in the transcriptomes of Plutella xylostella larvae infected 


and uninfected by saria fumosorosea 


GO 条 目 
GO term 

















差异 表达 基因 数量 
Number of DEUs 














细胞 组 分 Cellular component 
胞 外 区 域 Extracellular region 





分 子 功能 Molecular function 
氧化 还 原 酶 活性 Oxidoreductase activity 
细菌 细胞 表面 结合 Bacterial cell surface binding 
细胞 表面 结合 Cell surface binding 
大 链 内 切 酶 抑制 剂 活性 Endopeptidase inhibitor activity 
大 酶 抑制 剂 的 活性 Peptidase inhibitor activity 
大 链 内 切 酶 调节 活性 Endopeptidase regulator activity 












































氧化 还 原 酶 活性 ,作用 于 匹配 的 供 体 ,与 匹配 供 体 发 生 氧 化 反应 ,引起 氧 分 子 还 原 成 两 个 水 分 子 


Oxidoreductase activity, acting on paired donors, with oxidation of a pair of donors resulting in the reduction of molecular 2 








oxygen to two molecules of water 


肽 酶 调节 活性 Peptidase regulator activity 








生物 学 过 程 Biological process 


应 激 反应 Response to stress 


翻译 起 始 的 eIF2a 磷酸 化 调节 Regulation of translational initiation by eIF2 alpha phosphorylation 


防御 反应 Defense response 














E ARAKI Regulation of autophagy 
然 免疫 反应 调控 Regulation of innate immune response 











度 系 统 过 程 Immune system process 
然 免 疫 反 应 Innate immune response 
包 因 子 生成 调控 Regulation of cytokine production 

















E imcgk oum 


度 反 应 Immune response 


应 激 反应 中 对 翻译 的 调控 Regulation of translation in response to stress 
应 激 反应 中 对 翻译 起 始 的 调控 Regulation of translational initiation in response to stress 
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RI 健康 和 感染 玫 烟 色 棒 束 孢 的 小 菜 蛾 幼虫 转录 组 中 差异 表达 基因 (Q H <0.05) 的 KEGG pathway 富 集 分 析 
Table 3 Significantly enriched KEGG pathways of differentially expressed unigenes ( DEUs) (Q-value <0.05) in the 


transcriptomes of Plutella xylostella larvae infected and uninfected by Isaria fumosorosea 


























ig ID 通路 差异 表达 基因 数量 Qi 

Pathway ID Pathway Number of DEUs Q-value 
ko05164 流感 A Influenza A 20 3.31e - 06 
ko04974 蛋白 消化 和 吸收 Protein digestion and absorption 18 1.41e -05 
ko04972 胰腺 分 泌 Pancreatic secretion 21 1.50e -05 
ko00965 甜菜 红色 素 生 物 合成 Betalain biosynthesis 5 4. 62e -05 
ko00073 角质 、 软 木 脂 和 量 质 的 生物 合成 Cutin, suberine and wax biosynthesis 6 1.75e -04 
ko04146 过 氧化 物 酶 体 Peroxisome 12 8.22e - 04 
ko00950 异 唆 啉 类 生物 碱 的 生物 合成 Isoquinoline alkaloid biosynthesis 5 9.31e -04 
ko04080 神经 活性 的 配 体 - 受 体 互 作 Neuroactive ligand-receptor interaction 18 1.48e -03 
ko00740 核 黄 素 代谢 Riboflavin metabolism 4 2.277e -03 
ko04610 补体 和 凝血 级 联 反 应 Complement and coagulation cascades 6 7.21e -03 
ko05014 肌 蓉 缩 Amyotrophic lateral sclerosis ( ALS) 6 7.21e -03 
ko00630 乙 醛 酸 和 二 羚 酸 代谢 Glyoxylate and dicarboxylate metabolism 5 1.71e -02 
ko00350 酷 氨 酸 代 谢 Tyrosine metabolism 6 2.60e -02 


" RTT 4) V (pancreatic secretion) (13 个 上 调和 8 个 
3 DEUs) o 3x 3 个 通路 所 涉及 到 的 基因 分 别 占 
总 注释 基因 的 12.27% , 11. 0496 和 12. 8896 ,并 且 
上 调 表达 DEUs 的 数量 多 于 下 调 表达 的 数量 。 一 些 
DEUs 同时 参与 了 多 个 代谢 通路 , 比如 高 表达 上 调 
基因 CCG001601. 1 (log,Ratio 29. 1) ,编码 着 胰岛 
素 (trypsin) ,同时 参与 了 3 个 代谢 通路 ,与 虫 体内 能 
量 代谢 调控 密切 相关 。 此 外 , 富 集 DEUs 较 多 的 通 
路 还 有 "过 氧化 物 酶 体 ”(peroxisome) “神经 活性 的 
配 体 - 受 体 相互 作用 ”( neuroactive ligand-receptor 
interaction ) 等 ,涉及 传染 性 疾病 、 能 量 代谢 、 生 理 防 
御 和 神经 性 疾病 等 方面 。 


















































3 讨论 








昆虫 病原 真菌 是 昆虫 病原 微生物 中 最 大 的 一 个 
类 群 ,其 作为 害虫 生物 防治 的 重要 组 成 部 分 ,在 持续 
控制 害虫 .保护 生态 平衡 方面 具有 不 可 符 代 的 作用 。 
尽管 真菌 杀 虫 剂 安 全 无 环境 污染 ,但 与 化 学 农药 相 
比 , 却 存在 着 杀 虫 速率 缓慢 的 弱点 ( 王 滨 等 ,2001; 
Ren et al., 2011) 。 昆 虫 病原 真 戎 在 寄主 表皮 定 殖 
和 侵入 昆虫 血 腔 增殖 的 两 个 时 期 ,是 昆虫 免疫 识别 
与 防御 病原 真菌 对 抗 寄 主 免疫 反 应 的 关键 时 期 E 
虫 与 病原 真菌 之 间 识 别 与 反 识 别 、 免 疫 与 抗 免疫 的 
动态 互 作 关系 影响 着 真菌 的 作用 速率 。 

球 孢 白 僵 菌 和 金龟 子 绿 僵 苗 侵 染 昆虫 时 ,伴随 
着 白 僵 菌 素 (beauvericin ) 和 破坏 素 ( destruxin ) 等 毒 
KEPENK K RAEE, 2002) 。 目 前 , 玫 烟 色 

































































棒 束 孢 能 否 分 泌 毒 素 还 不 清楚 , REON Co ER T 
EBCL 03011 致 病 速 度 之 快 的 原因 仍 不 得 而 知 ( 吕 利 
华 等 , 2007) ,其 对 寄主 昆虫 的 致 病 机 制 有 可 能 与 依 
靠 毒 素 致死 宿主 的 病原 真菌 不 同 ,推测 小 菜 蛾 对 玫 
烟 色 棒 束 孢 的 免疫 方式 有 其 独特 性 。 本 研究 通过 
RNA-seq 技术 测定 了 小 菜 蛾 4 12:57] nt p BOR C pe 
ZR 12 h 后 的 表达 谱 , 与 健康 小 菜 蛾 表达 谱 数 据 
比 对 ,在 获得 的 差异 表达 基因 中 ,挖掘 了 一 批 与 玫 烟 
色 棒 束 孢 侵 染 直接 或 潜在 相关 的 小 菜 蛾 免疫 应 答 基 
因 与 相关 调控 网 路 。 在 感 贡 和 健康 小 菜 蛾 转录 组 这 
两 个 表达 谱 里 , 共 得 到 15 089 个 差异 表达 基因 ,其 
中 ,显著 差异 表达 351 个 ,上 调 表 达 基 因 275 个 ,下 
调 表达 基因 76 个 ,与 免疫 潜在 相关 的 基因 有 156 
个 ;GO 富 集 分 析 有 102 个 差异 基因 分 布 到 46 个 GO 
term 里 ,KEGG pathway 富 集 分 析 结 果 显 示 有 163 个 
差异 基因 分 布 到 153 个 代谢 通路 里 ,这 些 基因 中 , 包 
含 着 大 量 潜 在 的 与 免疫 调控 相关 的 基因 。 其 中 有 与 
病原 菌 识别 相关 的 1 个 B-1,3- 葡 聚 糖 识 别 和 蛋白 
( BGRP) ( CCG009703. 1) 1 个 模式 识别 受 体 - 革 兰 
氏 阴 性 菌 结合 蛋白 (GNBP) (CCG009706.1) 和 4 个 
受 体 - 肽 聚 糖 识别 蛋白 (PGRP) ( CC6015207. 1, 
CCG004942.1, CCGO008494. 1 和 CCG001311.2) 经 
BTA EIRA. PAMA Drosomycin 未 被 显 
著 诱 导 或 抑制 ,可 能 与 小 菜 蛾 对 真菌 免疫 特异 性 有 
关 , 也 可 能 与 取样 时 间 点 未 诱导 这 些 基 因 表 达 有 关 。 
已 有 报道 病原 真菌 可 诱导 Toll 信号 通路 ,经 模式 识 
别 受 体 识别 后 ,诱导 丝氨酸 蛋白 酶 级 联 反 应 ,进而 激 
活 Spaetzle 配 体 产生 构象 变化 以 结合 Toll. 受 体 
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( Bischoff et al., 2004) 。 本 研究 Toll 信号 途径 中 有 3 
个 Toll 受 体 (CCG003668. 1，CCG001789. 1 和 
CCG007264. 1) 被 诱导 显著 上 调 表达 ,Imd 途径 的 一 
些 基因 也 被 不 同 程度 地 诱导 表达 ,两 种 信号 通路 间 
并 非 相 互 孤立 (Kaneko and Silverman, 2005) , 感 菌 
虫 体 为 了 维持 机 体 的 生理 平衡 ,Imd 和 Toll 通路 上 
的 基因 被 激活 ,参与 抗菌 反应 。 

此 外 ,昆虫 经 病原 真菌 诱导 后 ,天然 免 疫 系 统 中 
还 会 出 现 类 似 抗 真菌 蛋白 的 分 泌 和 新 蛋白 的 产生 
(Wang et al., 2007) 。 昆 虫 的 储存 和 蛋白、 转 铁 蛋 白 、 
热 激 蛋 白 、 核 糖 体 相 关 膜 蛋白 .过 氧化 氢 酶 .NADH 
脱氧 酶 、 谷 胱 甘 肽 S- 转 移 酶 等 也 是 潜在 的 免疫 应 答 
相关 因子 ( Vierstraete et al., 2004; Hou et al., 
2011; 党 向 利 等 , 2012 ) 。 我 们 的 测序 结果 中 出 现 
了 一 些 参与 能 量 代 谢 、 核 酸 及 和 蛋白 代谢 和 加 工 、 胁 迫 
及 解毒 等 过 程 的 基因 ,其 作用 具有 网 络 性 的 特点 ,不 
同 基 因 之 间 通 过 相互 协同 或 相互 制约 ,构成 了 复杂 
的 免疫 调节 网 络 。 本 研究 获得 的 表达 谱 数 据 库 ,为 
探讨 小 菜 蛾 与 玫 烟 色 棒 束 孢 间 的 互 作 提 供 了 大 量 的 
免疫 候选 基因 ,关于 这 些 差 异 表达 基因 的 功能 和 参 
与 的 途径 及 其 在 小 菜 蛾 对 抗 玫 烟 色 棒 孢 侵 染 时 免疫 
应 答 作 用 ,还 需要 今后 进一步 深入 研究 。 
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